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£h : ABSTRACT 

ov 

CN , The EMBH model, previously developed using GRB 991216 as a prototype, is here applied to GRB 980425. 

We fit the luminosity observed in the 40-700 keV, 2-26 keV and 2-10 keV bands by the BeppoSAX satellite. 
J> ■ In addition we present a novel scenario in which the supernova SN1998bw is the outcome of an "induced 
gravitational collapse" triggered by GRB 980425, in agreement with the GRB-Supernova Time Sequence 
^ | (GSTS) paradigm (Ruffini et al. 2001c). A further outcome of this astrophysically exceptional sequence 
' of events is the formation of a young neutron star generated by the SN1998bw event. A coordinated 
. observational activity is recommended to further enlighten the underlying scenario of this most unique 
q I astrophysical system. 

O 1 INTRODUCTION 

The aim of this talk is to present the application of the EMBH theory, previously successfully applied to 
C3 '■ GRB 991216 used as a prototype, to the case of GRB 980425 (Pian et al, 2000) and SN1998bw (Galama et 
;> al., 1998). This is a particularly important test for the validity of the EMBH theory over a range of energies 
of 6 orders of magnitude: as we will see, both sources appear to be spherically symmetric and the respective 
total energies are Etot — 5 x 10 53 ergs and Etot — 10 48 ergs. We recall that the EMBH theory (see Ruffini 
et al., 2003b) depends only on three parameters, the energy of the dyadosphere Edy a , the B parameter 
and the factor 1Z describing the interstellar medium (ISM) porosity. The theory, therefore, explains all the 
observed features of the bolometric intensity variations of the afterglow as well as the spectral properties 
of the source and, in the specific case of GRB 980425, it also allows to clarify the general astrophysical 
scenario in which the GRB actually occurs. In this system, in fact, we propose that GRB 980425 has 
been the trigger of a phenomenon of "induced gravitational collapse" (Ruffini et al., 2001c) originating the 
supernova explosion and we also witness the birth of a young neutron star out of the supernova event. 
This extraordinary coincidence of these three astrophysical events represents an unprecedented scenario of 
fundamental importance in the field of relativistic astrophysics. Using the EMBH theory we shall explore: a) 
the process of black hole formation in the event GRB 980425 (Pian et al., 2000), b) the concept of "induced 
gravitational collapse", introduced in the GRB-Supernova Time Sequence (GSTS) paradigm (Ruffini et al., 
2001c), and its link to a very special supernova type in SN1998bw (Galama et al., 1998), and finally c) the 
observation for the first time of the cooling of a hot newly formed neutron star. 

The observational situation of this system is quite complex. In addition to the source GRB 980425 
and the supernova SN1998bw, two X-ray sources have been found by BeppoSAX in the error box for the 
location of GRB 980425: a source SI and a source S2 (Pian et al., 2000). Since the nature of the two 
sources SI and S2 was not clear, a variety of slopes in the decaying part of the afterglow have been proposed 
(see Fig. 1). Kulkarni et al. (1998) have proposed to explain both the supernova SN1998bw and the 
GRB 980425 observations by a new class of GRBs, distinctly different from the cosmological ones, both 
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originated by a single unusual supernova event. Similarly, Iwamoto et al. (1998) have tried to explain both 
the supernova event and the GRB with a new kind of supernova with an extremely large explosion energy, 
larger than 10 52 ergs, which they identify with the "hypernovae" predicted by Paczyhski (1998). In this 
approach a totally novel concept is introduced: the supernova itself is assumed to originate in the process 
of gravitational collapse to a black hole of a massive progenitor star (~ 4OM0) with a particularly large 
angular momentum and strong magnetic field. A large rotational energy of the black hole extracted with a 
strong magnetic field is called in, by these authors, to explain the successful explosion of this "hypernova" 
leading both to the GRB and the supernova. 

Our approach is drastically different. We first interpret the GRB 980425 within the EMBH theory. 
This allows the computation of the luminosity, spectra, Lorentz gamma factors, and more generally all 
the dynamical aspects of the source. Having characterized the features of GRB 980425, we can gradually 
approach the remaining part of the scenario, disentangling the GRB observations from the supernova ones 
and from the sources SI and S2. This leads to a natural time sequence of events and to their autonomous 
astrophysical characterization. 

THE ENERGETICS, DYNAMICAL PARAMETERS AND SPACE-TIME PARAMETRIZA- 
TION OF GRB 980425 

Our approach has focused on identifying the energy extraction process from the black hole (Christodoulou 
& Ruffini, 1971) as the basic energy source for the GRB phenomenon. The distinguishing feature is a 
theoretically predicted source energetics all the way up to 1.8 x 10 54 (Mbh/Mq) ergs for 3.2M < M BH < 
7.2 x 10 6 Mq (Damour &; Ruffini, 1975). In particular, the formation of a "dyadosphere" , during the 
gravitational collapse leading to a black hole endowed with electromagnetic structure (EMBH) has been 
indicated as the initial boundary conditions of the GRB process (Ruffini, 1998; Preparata et al., 1998). Our 
model has been referred as "the EMBH model", although the EMBH physics only determines the initial 
boundary conditions of the GRB process by specifying the physical parameters and spatial extension of the 
neutral electron positron plasma originating the phenomenon created in the dyadosphere. The creation of 
this plasma is due to the vacuum polarization process occurring in a supercritical field by the Heisenberg- 
Euler-Schwinger process (see Heisenberg & Euler, 1935; Schwinger, 1951; Damour & Ruffini, 1975; Preparata 
et al, 1998). 

Traditionally, following the observations of the Vela (Strong, 1975) and CGRO 1 satellites, GRBs have 
been characterized by few parameters such as the fluence, the characteristic duration (T90 or T50) and the 
global time averaged spectral distribution (Band et al., 1993). With the observations of BeppoSAX 2 and 
the discovery of the afterglow, and the consequent optical identification, the distance of the GRB source has 
been determined and consequently the total energetics of the source has been added as a crucial parameter. 

The observed energetics of GRBs, computed for spherically symmetric explosions, do coincide with the 
ones theoretically predicted in Damour & Ruffini (1975). This has been the major reason which has moti- 
vated us to reconsider and develop in full details the EMBH model. For simplicity, we have considered the 
vacuum polarization process occurring in an already formed Reissner-Nordstrom black hole (Ruffini, 1998; 
Preparata et al., 1998), whose dyadosphere has an energy Ed ya . It is clear, however, that this is only an 
approximation to the real dynamical description of the process of gravitational collapse to an EMBH. In 
order to prepare the background for attacking this extremely complex dynamical process, we have clari- 
fied some basic theoretical issues, necessary to implement the description of the fully dynamical process of 
gravitational collapse to an EMBH (see Ruffini & Vitagliano, 2002,2003; Cherubini et al., 2002). 

We have then given the constitutive equations for the five eras in the EMBH model (see for details Ruffini 
et al., 2003a and references therein). The Era I: the e + e~ pairs plasma, initially at 7 = 1, self propels itself 
away from the dyadosphere as a sharp pulse (the PEM pulse), reaching Lorentz gamma factor of the order 
of 100 (Ruffini et al., 1999). The Era II: the PEM pulse, still optically thick, engulfs the remnant left over in 
the process of gravitational collapse of the progenitor star with a drastic reduction of the gamma factor; the 
mass Mb of this engulfed baryonic material is expressed by the dimensionless parameter B = Mbc 2 / Ed ya 
(Ruffini et al., 2000). The Era III: the newly formed pair-electromagnetic-baryonic (PEMB) pulse, composed 



http://cossc.gsfc.nasa.gov/batse/ 
http: / /www.asdc.asi.it /bepposax/| 



3 



10 



Detector arrival time (t a ) (years) 



10 



-8 



10 



-7 



io- 



10" 



10" 



10 



-3 



10 



-2 



10 



1 



44 



42 



CO 

i 5 10 40 



WFC data 2-10keV 
WFC u.l. 2-10keV 
MECS data 2-10keV for S1 
MECS data 2-10keV for S2 

-1.6 
j-1.4 



10 L 



10 



10" 



-7 



- 10" 



- 10" 



- 10 



38 



10" 11 ~° 

' u CD 

> 



- 10" 



10 



36 



- 10" 



10 



34 



10" 



10" 



10' 



10' 



10" 



10 H 



10~ 



10 c 



10' 



Detector arrival time (t a ) (s) 



10 c 



Fig. 1. The BeppoSAX MECS light curves in the 2-10 keV band of the X-ray sources SI and S2 from Pian et al. 
(2000) detected in the GRB 980425 fields are shown. The WFC early measurement and 3a upper limits in the 
same band are also shown. The dotted lines represent the power laws of indexes ~ 1.4 and ~ 1.6 connecting the 
last WFC measurement and the first and last of the first 5 MECS points of S2 which are in this picture replaced 
by the first point that has been obtained by integrating and averaging the flux in shorter time interval, as in Pian 
et al. (2000). In the vertical axis is reported (left) the luminosity of the source and (right) the observed flux as 
a function of the arrival time at the detector in seconds (bottom) and in years (top). 
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of e + e~ pair and of the electrons and baryons of the engulfed material, self-propels itself outward reaching in 
some sources Lorentz gamma factors of 10 3 -10 4 ; this era stops when the transparency condition is reached 
and the emission of the proper-GRB (P-GRB) occurs (Bianco et al., 2001). The Era IV: the resulting 
accelerated baryonic matter (ABM) pulse, ballistically expanding after the transparency condition has been 
reached, collides at ultrarelativistic velocities with the baryons and electrons of the interstellar matter (ISM) 
which is assumed to have an average constant number density, giving origin to the afterglow. The Era V: 
this era represents the transition from the ultrarelativistic regime to the relativistic and then to the non 
relativistic ones (see Ruffini et al., 2003a, 2003b and references therein). 

The EMBH model differs in many respects from the models in the current literature. The major difference 
consists in the following points: 

a) The appropriate theoretical description of all the above mentioned five eras is implemented, as well as 
the evaluation of the process of vacuum polarization originating the dyadosphere. The description of the 
inner engine originating the GRBs has never been addressed in the necessary details in the literature; 

b) The dynamical equations as well as the description of the phenomenon in the laboratory time and the time 
sequence carried by light signals recorded at the detector have been explicitly integrated (see e.g. Ruffini 
et al., 2003a, 2003b). In doing so we have also corrected a basic conceptual inadequacy, common to all 
the current works on GRBs, which led to an improper spacetime parametrization of the GRB phenomenon, 
preempting all these works from their predictive power: the relation between the photon arrival time at 
the detector and their emission time in the laboratory frame, expressed in our approach by an integral 
of a function of the Lorentz gamma factor extended over all the past source worldlines, has been in the 
current literature expressed as a function of an instantaneous value of the Lorentz gamma factor. These 
two approaches are conceptually very different and lead to significant qualitative differences (Ruffini et al., 
2003a, 2003b and references therein). 

c) The treatment of the afterglow has been widely considered in the literature by the so-called "fireball 
model" (see e.g. Meszaros & Rees, 1992, 1993; Rees & Meszaros, 1994; Piran, 1999 and references therein). 
However, also in the description of the afterglow, there are major differences between the works in the 
literature and our approach (Ruffini et al., 2003b). 

The equations of motion in our model depend only on two free parameters: the total energy Etot, which 
coincides with the dyadosphere energy E dya , and the amount Mb of baryonic matter left over from the 
gravitational collapse of the progenitor star, which is determined by the dimensionless parameter B = 
M B c 2 /E dya . The best fit of GRB 980425 is reproduced in Tab. 1. It correspond to E dya = 1.1 x 10 48 ergs, 
B = 7xl0~ 3 and the ISM average density is found to be (rii sm ) = 0.02particle/cm 3 . The plasma temperature 
and the total number of pairs in the dyadosphere are respectively T = 1.028 MeV and N e ± = 5.3274 x 10 53 . 

THE GRB 980425 LUMINOSITY IN SELECTED ENERGY BANDS PREDICTED BY THE 
EMBH MODEL 

Recently, within the EMBH model, we have developed an attempt to theoretically derive the GRB spectra 
out of first principles as well as the GRB luminosity in fixed energy bands (Ruffini et al. 2003c) . We have 
adopted three basic assumptions: a) the resulting radiation as viewed in the comoving frame during the 
afterglow phase has a thermal spectrum and b) the ISM swept up by the front of the shock wave, with a 
Lorentz gamma factor between 300 and 2, is responsible for this thermal emission, c) We also adopt, like 
in our previous papers (Ruffini et al., 2001a,2001b, 2002, 2003b), that the expansion occurs with spherical 
symmetry. This three assumptions are different from the ones adopted in the GRB literature, in which the 
afterglow emission is believed to originate from synchrotron emission in the production of a shock or reverse 
shock generated when the assumed jet-like ejecta encounter the external medium (see e.g. Giblin et al., 
2002 and references therein). 

In the EMBH model the structure of the shock is determined by mass, momentum and energy conservation: 
the constancy of the specific enthalpy, which is a standard condition in shock rest frames (Zel'dovich & 
Rayzer, 1966) and have been used in our derivation (Ruffini et al., 2003b). The only free parameter of our 
model is the size of the "effective emitting area" in the shock wave front: A e ff. Since the determination 
of this free parameter is performed here by empirically fitting the observational data, we avoid ambiguities 
due to the absence of relevant theoretical and laboratory results on relativistic shocks for Lorentz factor 
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Table 1. Gamma factor for selected events and their space-time coordinates: r is the radial coordinate in the 
laboratory frame, t is the laboratory time and t d a is the photon arrival time at the detector. The points 1,2,3,4 
mark the beginning of first four "eras" (see Ruffini et al., 2003a, 2003b) until the transparency point (point 4). 
The points A, B,C correspond to three events respectively before, during and after the peak of the E-APE. 
The points D, E, F,G, H, I, J are taken in the decaying part of the afterglow. The last column shows how the 
apparent motion in the radial coordinate, evaluated in the arrival time at the detector, leads to an enormous 
apparent "superluminal" behaviour. This illustrates the impossibility of using such a classical estimate in regimes 
of high gamma Lorentz factor. 
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7 ~ 300. 

The temperature T of the black body in the comoving frame is then 



( A ^ V /4 (i) 



where 



K = ^ (2) 

is the ratio between the "effective emitting area" and the ABM pulse surface A a i> m , a is the Stefan-Boltzmann 
constant and AE[ nt is the proper internal energy developed in the collision between the ABM pulse and the 
ISM in a proper time interval At (see Ruffini et al., 2003b, 2003c). The ratio 7Z, which is a priori a function 
that varies as the system evolves, is evaluated at every given value of the laboratory time t. 

All the subsequent steps are now uniquely determined by the equations of motion of the system. The 
basic tool in this calculation involves the definition of the EQuiTemporal Surfaces (EQTS) for the relativistic 
expanding ABM pulse as seen by an asymptotic observer. The key to determining such EQTS (see Fig. 1 in 
Ruffini et al., 2002) is the relation between the time t in the laboratory frame at which a photon is emitted 
from the ABM pulse external surface and the arrival time at which it reaches the detector. We have 
instead adopted the equations (see Ruffini et al., 2002,2003b): 

tj. (i + ,) (t - r M cos „ + _ (1 + z) L _ /.' » C)f + *■ „, „ + g j 

= (1 + z) (t-COS#£ ^1 - ^2j^j dt ' + ^ (1 - COS §) \ , (3) 
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Fig. 2. a) The light curve computed within the EMBH model in 7-ray (40-700 keV) is represented; the 
horizontal bar represents the peak flux in the 40-700 keV band measured by the GRBM (Frontera et al., 2000). 
The horizontal dotted line represents the noise level of the GRBM detector, b) The light curve computed within 
the EMBH model in hard X-rays in the 2-26 keV band with the peak flux and time duration from WFC (Frontera 
et al., 2000). c) The light curve computed within the EMBH model in hard X-rays in the 2-10 keV band with 
the experimental data in the same band from WFC (Pian et al., 2000) 

where z is the cosmological redshift of the source, ro = r (t = 0) and $ is the angle subtended by the emission 
point of the photon on the ABM pulse external surface, having defined $ = along the line of sight. For such 
a relation, different approximations exist in the literature (see e.g. Fenimore et al., 1996; Sari, 1997,1998, 
Rees & Meszaros, 1998; Fenimore et al., 1999; Granot et al., 1999 and see also the reviews by Piran, 1999; 
Meszaros, 2002). There are both quantitative and qualitative important differences in the use of Eq.(3) 
instead of the ones usually adopted in the literature (Ruffini et al., 2003b). It is a matter of fact that only 
the use of the EQTS defined by Eq.(3) allows to fit the observational data and that even a very minor 
departure from it leads to unacceptable results. It is also important to recall that this inadequacy in the 
current literature in the relation between the arrival time and the laboratory time has affected also all the 
estimates of the power-law slopes in the afterglow, preempting all current theoretical considerations in the 
literature of their predictive power. All the considerations about beaming in GRBs existing in the current 
literature have to be reformulated on the ground of the proper theoretical treatment and of Eq.(3) (see 
Ruffini et al., 2003b). Having so determined the EQTS we have been able to evaluate the source luminosity 
in a given energy band, in agreement with the above mentioned new assumptions (Ruffini et al., 2003c). 

We can now proceed to the best fit of the GRB 980425 observed data. The best fit of the observed 
luminosity in selected energy bands has been obtained for the above mentioned values of the parameters 
Edya, B and for 1Z monotonically decreasing from 4.81 x 10~ 10 to 2.65 x 10~ 12 . The results are given in 
Fig. 2 where the luminosity is computed as a function of the arrival time for three selected energy bands. 

We can then conclude: 

1. The best fit is obtainable under almost perfect spherical symmetry. This has been proven as a result of 
an analysis of unprecedented redundancy: the luminosity curves are obtained from an integration over 
almost 10 7 different paths, relating the observer to the EQTS. This procedure tests, to a very high 
level of accuracy, any departure from spherical symmetry as well as any departure from the computed 
equations of motion of the source (Ruffini et al., 2003b). This same circumstance was encountered in 
GRB 991216 (Ruffini et al., 2003b). 

2. Each luminosity curve as a function of the arrival time presents complex behavior, which could be 
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Fig. 3. The light curves in selected bands are reported as well as the MECS light curves in the 2-10 keV band 
of SI and S2 (Pian et al. 2000) as well as the optical data (Iwamoto, 1999). Here are also reported theoretical 
models of neutron star cooling (Canuto, 1978). 



erroneously interpreted as evidence for breaks in the power-law indexes leading to erroneous inferences 
on the possible existence of jets (Ruffini et al., 2003c). 

3. Although each luminosity curve presents some special features, the bolometric luminosity has a very 
clear and simple power-law behavior with the "golden value" index n = — 1.6 (Ruffini et al., 2003b). 

THE NEW ASTROPHYSICAL SCENARIO AND THE NEWLY BORN NEUTRON STAR 

In Fig. 3 the luminosities in the three bands are represented together with the optical data of SN1998bw 
(black dots), the source SI (black squares) and the source S2 (open circles). It is then clear that GRB 980425 
is separated both from the supernova data and from the sources SI and S2. 

While the occurrence of the supernova in relation to the GRB has already been discussed with the 
GRB-Supernova Time Sequence (GSTS) paradigm (Ruffini et al. 2001c), we like to address here a different 
fundamental issue: the possibility of observing the birth of a newly formed neutron star out of the supernova 
event, which in turn has been triggered by the GRB 980425. 

In the early days of neutron star physics it was clearly shown by Gamow & Schoenberg (1940,1941) 
that the URCA processes are at the very heart of the supernova explosions. The neutrino-antineutrino 
emission described in the URCA process is the essential cooling mechanism necessary for the occurrence of 
the process of gravitational collapse of the imploding core. Since then, it has become clear that the newly 
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formed neutron star can be still significantly hot and in its early stages will be associated to three major 
radiating processes (Tsuruta, 1964,1979,2002; Canuto, 1978): a) the thermal radiation from the surface, b) 
the radiation due to neutrino, kaon, pion cooling, and c) the possible influence in both these processes of 
the superfluid nature of the supra-nuclear density neutron gas. Qualitative representative curves for these 
cooling processes, which are still today very undetermined due to the lack of observational data, are shown 
in Fig. 3. 

It is of paramount importance to follow the further time history of the two sources SI and S2. If, as 
we propose, S2 is a background source, its flux should be practically constant in time and this source has 
nothing to do with the GRB 980425 / SN1998bw system. If SI is indeed the cooling radiation emitted by 
the newly born neutron star, it should be possible to notice a very drastic behavior in its luminosity as 
qualitatively expresses in Fig. 3. 

CONCLUSIONS 

It is particularly attractive, in conclusion, to emphasize some of the analogies and the differences between 
the case of GRB 991216 and the one of GRB 980425: 

1. In both these sources a GRB and, independently, a supernova event are present. In the case of 
GRB 991216 the inferences of the supernova can be obtained only on the ground of the emission 
of iron lines (Piro et al., 2000; Ruffini et al., 2001c). In the present case of GRB 980425 we have 
a very fortunate circumstance: the GRB source is much weaker and is much closer to us (z = 1.0 
for GRB 991216 and z = 0.00835 for GRB 980425). This situation is particularly important for 
obtaining detailed data on the supernova and on the possible occurrence of a newly born neutron 
star. This occurrence could not be observed in GRB 991216 due to the very large distance and to the 
overwhelming X-ray luminosity of the afterglow (see Ruffini et al. , 2003b) . 

2. The energetics of GRB 991216 is E dya = 4.83 x 10 53 ergs, while the one of GRB 980425 is E dya = 
1.1 x 10 48 ergs. It is very impressive that the EMBH model applies over a range of more than 6 
orders of magnitude, giving important inferences on the source as well as on the structure of the 
ISM surrounding the source. It is significant that in both sources the condition of spherical symmetry 
appears to be strongly implemented and this fact is a very clear discriminant among all possible sources 
of energy for GRBs. 

3. The main difference between GRB 991216 and GRB 980425 is then traced back, within the EMBH 
model, to the different parameters occurring in the dyadosphere and to the nature of the "effective" 
Reissner-Nordstrom geometry which we have used as a reliable estimate for the dynamical processes of 
gravitational collapse leading to the formation of the EMBH. For GRB 991216 we have a ratio between 
the electromagnetic energy and the total mass of the imploding core, described by the parameter Q/M, 
given by Q/M = 0.23 while for GRB 980425 we have Q/M = 6.5 x 10~ 4 . In both cases a reasonable 
mass of the black hole is M BH ~ 10M Q (Ruffini et al., 2003b). 

A dedicated observational campaign, both with XMM and Chandra, to follow the cooling of the newly 
formed neutron star is needed in order to gain for the first time information on this extremely important 
astrophysical process. 
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